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There is a growing consensus that physics majors need to learn computational skills, but many 
departments are still devoid of computation in their physics curriculum. Some departments may 
lack the resources or commitment to create a dedicated course or program in computational physics. 
One way around this difficulty is to include computation in a standard upper-level physics course. 
An intermediate classical mechanics course is particularly well suited for including computation. 
We discuss the ways we have used computation in our classical mechanics courses, focusing on how 
computational work can improve students' understanding of physics as well as their computational 
skills. We present examples of computational problems that serve these two purposes. In addition, 
we provide information about resources for instructors who would like to include computation in 
their courses. 



PACS numbers: D01.30.Pp, S01.40.Fk, N01.50.H- 
I. INTRODUCTION 

The primary purpose of this article is to suggest a 
method for incorporating computation into upper-level 
classical mechanics courses. There is an emerging consen- 
sus in the physics community that computational skills 
are important for physicists,^ but too often there is little 
or no computation included in the physics curriculumP 
Probably the best way to strengthen the computational 
component of the physics curriculum is to incorporate 
computation into all components of the curriculum. Sev- 
eral institutions have developed large-scale computa- 
tional physics programs,^ which typically include at least 
one dedicated computational physics course as well as 
computational components in other upper-level courses. 
Other institutions offer a single computational physics 
course along with an otherwise traditional curriculum. 
With a variety of good computational physics textbooks 4 
available it may appear that there is no reason not to offer 
at least one course in computational physics. Unfortu- 
nately, the reality is that there are constraints that may 
prevent some departments from offering such a course. 
These constraints may involve the number of physics 
faculty, the computational background of the physics 
faculty, low student enrollment in physics courses, or 
even political resistance to curricular changes. For fac- 
ulty who face these constraints but desire to include 
some computation in the physics curriculum, the best 
approach may be to incorporate computation into exist- 
ing courses P Even departments that have a computa- 
tional physics course may be looking for ways to increase 
the role of computation in the physics curriculum. We 
believe that the intermediate classical mechanics course, 
typically taken by students in their sophomore or junior 
year, is ideally suited for incorporating computation. 

Computation can contribute to a variety of learning 
goals, and different approaches to incorporating compu- 



tation into physics courses will likely emphasize differ- 
ent goals. One goal of including computation in physics 
courses is simply to improve students' learning of physics 
concepts. Computer visualizations and interactive simu- 
lations are particularly useful in this regard. Computa- 
tion can also open the door to new and important topics 
such as nonlinear dynamics and chaos. Including com- 
putation in physics courses can also increase physics stu- 
dents' familiarity with widely-used computational tools 
and help students see how these tools can be applied to 
solving physics problems. Finally, computation can be 
used to introduce students to important numerical al- 
gorithms. In many cases these algorithms can provide 
insight into important physics concepts, in addition to 
serving as tools for carrying out computations. 

Although all of these learning goals can be achieved 
through a dedicated course in computational physics, 
as mentioned above this may not be a practical ap- 
proach in all departments. Incorporating computation 
into a standard course may be an effective way to in- 
troduce computation into the physics curriculum, per- 
haps as a temporary solution until a dedicated compu- 
tational course can be added. Computation can be in- 
cluded in a standard course even if the students have 
no computational background (although, of course, more 
can be done if the students have had a course in pro- 
gramming or computational physics). We feel that the 
standard (sophomore/junior level) intermediate classical 
mechanics course is well suited for including computa- 
tion. Students in classical mechanics can benefit tremen- 
dously from computer visualizations, which help them to 
build intuition about classical dynamics. Classical me- 
chanics provides an excellent forum for introducing a va- 
riety of important numerical tools such as ODE solvers, 
root finding, numerical integration, numerical linear al- 
gebra, etc. Some important topics in modern classical 
mechanics, like chaos, cannot be effectively taught with- 
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out computation. Another advantage of the classical me- 
chanics course is that it is typically taught early in the 
upper-level physics curriculum, often immediately after 
the introductory sequence. Introducing students to com- 
putation at this early stage gives them the opportunity to 
use their computational skills in later physics and math- 
ematics courses. 



II. IMPLEMENTATION ISSUES 

To begin using computation in a classical mechanics 
course one must first make choices about which com- 
putational tools to use and how to use them. Choos- 
ing the right computational platform can be difficult, as 
each platform has advantages and disadvantages. Sym- 
bolic and numeric mathematics software packages (such 
as Mathematica, Maple, and Matlab, each of which has 
its own advantages and disadvantages®) offer quick and 
relatively easy ways to perform computational work, pro- 
vide a wide range of computational tools, and include 
high-quality visualization tools. The downside is that 
these tools can be used with little understanding of the 
numerical methods that are being employed. Further- 
more, these software packages can be quite expensive. 
The other main alternative is to have students create 
programs from scratch using a standard programming 
language like Fortran, Java, or Python. Students who 
do their own coding must develop an understanding of 
the models they are studying as well as the numerical 
methods that are employed. This solution is generally 
inexpensive since free compilers are available for these 
languages on a variety of operating systems. Many com- 
mon programming tasks can be carried out with the aid 
of freely-available numeric libraries, like Open Source 
Physics (OSP)Pand SciPyPto reduce the time that must 
be spent creating the programs. However, in spite of 
these libraries there is still a great deal of programming 
"overhead" that must be addressed and instructors whose 
students have no programming background may not have 
sufficient time to teach the necessary programming skills. 
There are some intermediate solutions available, such as 
the Easy Java Simulations (EJS) package, 9 which auto- 
mates many programming tasks and allows users to fo- 
cus on the numerical algorithms that are used. However, 
EJS does not offer the broad range of functionality that 
is available with a standard programming language or a 
software package like Mathematica. 

We have chosen to use Mathematica and Matlab in 
our classical mechanics courses, but we plan to move to- 
ward using Java/OSP and EJS for some topics. We like 
Mathematica for its powerful symbolic mathematics ca- 
pabilities, and Matlab because it is widely used in in- 
dustry. Both of these packages are tools that students 
will likely be able to use in other classes and perhaps 
throughout their careers. Programming in Java with the 
OSP library provides an opportunity to share both pro- 
grams and curricular materials with the broader physics 



community and contribute to (as well as benefit from) 
a larger project.^ EJS provides simple tools for quickly 
creating visualizations and simulations. Each of these op- 
tions, as well as the others listed above, will be more or 
less suited to an individual instructor's specific situation. 
Instructors must be guided by the background of their 
students (Do they have programming experience?), their 
own computational experience (With what platforms is 
the instructor familiar?), as well as other practical con- 
siderations (Are there funds available to pay for software 
licenses? Will the software be supported by qualified 
campus staff?). 

The choice of platform (or platforms) should also be 
guided by how the instructor intends to use the com- 
putational tools in the course. In our classical mechan- 
ics courses we use computation in two main ways: for 
in-class demonstrations and for student computational 
projects. Typically we will demonstrate computer so- 
lutions to physics problems as part of our in-class lec- 
ture. The code used to construct the solution is made 
available to students so that they can examine the code 
and "experiment" with it at their leisure. We then fol- 
low up these demonstrations by assigning computational 
projects that students must complete on their own or in 
small groups. Early in the course these projects may in- 
volve only minor modifications of the code used in the 
demonstrations, perhaps to study the same phenomenon 
in a new system. In this way the demonstration code 
serves as a template that helps the students complete 
the computational project. As the course proceeds we 
demand more from our students, expecting them to con- 
struct computational solutions without a template (but 
using computational tools they have seen before) . These 
computational projects can easily be made the basis 
of formal writing assignments in which students must 
present their analytic and numerical work along with 
figures, typeset equations (possibly in /LaTcX), and a 
thorough discussion of the physics. 



III. ALGORITHMS AND PHYSICS 

Because the software packages we use tend to hide the 
details of the numerical algorithms they employ, we feel 
that it is important to provide students with some ex- 
plicit instruction on algorithms. We focus on simple al- 
gorithms, rather than the sophisticated algorithms em- 
ployed by the software packages we use, for a few im- 
portant computational tasks. This instruction enhances 
students' knowledge of computation because it empha- 
sizes the fact that computations require some algorithms 
and that the choice of algorithm can critically affect the 
success of the computation. In addition, teaching stu- 
dents about algorithms can sometimes lead to a better 
understanding of important physics concepts. 

As an example of how we use algorithms in our course, 
let us consider an object of mass m oscillating in one di- 
mension on an ideal spring with force constant k. The 
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FIG. 1: Comparison of (a) Euler and (b) Euler-Cromer al- 
gorithms applied to the simple harmonic oscillator. In both 
cases k = 1 N/m, m = 1 kg, x(0) = 1 m, and v(0) = 0. 
The dashed line shows the exact orbit in phase space. The 
solid line shows the solution generated by the algorithm. Also 
shown is a cluster of trajectories initially distributed at ran- 
dom in a square centered on x = 1 m, v = 0. The locations 
of these trajectories resulting from each algorithm are shown 
at t = 0, 1.8, 3.6, and 5.4 s. 



equation of motion for this simple harmonic oscillator is 
easy to solve analytically, but we can also take a numer- 
ical approach to the problem. The simplest algorithm 
that we might consider for this purpose is the Euler al- 
gorithm: 



x n + v n At, 
v n + F n At/m, 



(1) 



where x n , v n , and F n represent the displacement from 
equilibrium, the velocity, and the net force at time t n = 
to + nAt. For this system F n — —kx n . Although the 
software packages we use have their own built-in ODE 
solvers, we can also use these packages to implement a 
simple algorithm like the Euler algorithm. The solid red 
curve in Figure [TJ a) shows the trajectory of the object in 
phase space (position versus velocity) generated by the 
Euler algorithm when m — 1 kg, k = 1 N/m, x(0) = 1 m, 
v(0) = 0, and At = 0.2 s. The dashed blue curve in Fig. 
[Ija) shows the exact solution. It is clear from this result 
that the Euler algorithm is unstable since the trajectory 
produced by the algorithm spirals continually outward 
away from the exact solution. This also means that the 
Euler algorithm does not conserve energy: it is clear that 
with the Euler algorithm the object's maximum displace- 
ment and maximum speed are steadily increasing, so its 
total mechanical energy must also be steadily increasing. 

The failure of the Euler algorithm to conserve energy 
is related to yet another failure: the failure to preserve 
phase space volume. According to Liouville's Theorem, 
a conservative system (such as our simple harmonic os- 
cillator) must preserve the volume (or area) occupied by 
an ensemble of trajectories in phase space. Figure [ija) 
shows the locations generated by the Euler algorithm for 
a trajectory ensemble at times t = 0, 1.8, 3.6, and 5.4 
s. At t = the points are randomly distributed within 
a square region around x = 1 m, v = 0. As time passes 



the cluster of points moves through phase space and the 
shape of the region becomes distorted. However, it is 
clear from Figure [TJa) that the size of the region also 
grows over time, in clear contradiction to Liouville's The- 
orem. In fact, it is easy to show that the Euler algorithm 
doesn't preserve phase space volume. If we treat Equa- 
tion [T] as a two-dimensional map, then this map will pre- 
serve phase space volume if and only if the Jacobian of 
the map, defined by 



J = 



0x n 



(2) 



has determinant equal to one. A quick calculation will 
show that | J | 7^ 1 for the Euler algorithm, in the case of 
the simple harmonic oscillator. 

These problems with the Euler algorithm can be fixed 
with a simple modification, leading to what is known 
as the Euler-Cromer algorithm. 11 This algorithm simply 
updates the velocity first, and then uses the new velocity 
to update the position. The equations describing this 
algorithm are: 



x i+ i — Xi + v i+1 At 
v l+1 — Vi + FiAt/m, 



(3) 



where all quantities are defined as for the Euler algo- 
rithm above. This algorithm is stable and conserves en- 
ergy, on average, for oscillatory motion!^ We can quickly 
illustrate this by applying this new algorithm to our sim- 
ple harmonic oscillator. The results are shown in Figure 
[ljb). It is clear that the solution produced by the algo- 
rithm remains close to the exact solution, coinciding with 
the exact solution at every quarter period. Similarly, the 
total energy oscillates about the correct value with a pe- 
riod equal to half of the oscillator's period. Furthermore, 
this algorithm appears to preserve phase space volume. 
This can be easily proved by showing that the Jacobian 
for Equation [3] has determinant equal to one. 

Once these algorithms (and their differences) have 
been demonstrated, students should be given a chance 
to use them. They can explore the behavior of the algo- 
rithms as the time step is increased or decreased. They 
can compare the results of these algorithms with the re- 
sults obtained using the built-in ODE solver supplied by 
the software package. They can make modifications to 
the model by adding drag forces or driving forces. Note 
that students are not restricted to drag forces linear in 
the velocity, or sinusoidal driving forces. In fact, it is very 
instructive to consider the case of a quadratic drag force 
(a nonlinear system) and compare the results to those ob- 
tained with a linear drag force (a linear system). If these 
non-conservative forces are added then neither algorithm 
will preserve phase-space volume, as can be easily seen 
by finding the determinant of the Jacobian (and noting 
that F n is no longer a function of x n only) . Carrying out 
this calculation using the Euler-Cromer algorithm makes 
it clear that for drag forces, where dF/dv < 0, the deter- 
minant of the Jacobian is less than one indicating that 
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phase space volume will shrink over time. This is exactly 
what one should expect for a dissipative system. 

Other algorithms for common numerical tasks may 
be worth discussing in class. Algorithms for root find- 
ing (Newton-Raphson method, bisection method) and 
numerical integration (trapezoid approximation, Simp- 
son's Rule, Monte Carlo methods) are useful and simple 
enough to present without absorbing too much class time. 
Although class time spent discussing algorithms may take 
away from time spent discussing physics, the above ex- 
ample shows that it is possible to teach important physi- 
cal concepts along with algorithms. As another example, 
the Newton-Raphson method illustrates the concept of 
a stable attractor which is important for understanding 
the dynamics of dissipative systems. 



IV. COMPUTATIONAL PROJECTS 




Computational assignments should be chosen with care 
to ensure that they take full advantage of what comput- 
ing offers. One way to do this is to use computation to 
break out of the restrictions imposed by the need for an 
analytic solution (or approximate solution). We have al- 
ready mentioned the case of the harmonic oscillator with 
a quadratic drag force. Two other examples of this type 
are motion in a non-inertial reference frame and motion 
of a charged particle in electric and magnetic fields. Tra- 
ditional courses tend to focus on cases where analytic 
solutions are possible, such as the Eastward deflection of 
an object that drops from rest to Earth or the cyclotron 
motion of an electron in a magnetic field. Computation 
allows students to solve more realistic problems such as 
long-range projectile motion on Earth and the motion 
of an electron in combined electric and magnetic fields. 
Solving such problems computationally not only allows 
students to tackle realistic problems, it also allows stu- 
dents to visualize the motion in these systems. Figure [2] 
shows two calculated trajectories, one with non-inertial 
forces and one without, for a long-range projectile fired 
due East from Rome, Georgia. Visualizing this motion 
helps students get a sense of the magnitude of the effect 
of non-inertial forces. Figure [3] shows an OSP application 
that simulates a charged particle moving in the presence 
of constant (but general) three dimensional electric and 
magnetic fields. Students can experiment with the pa- 
rameters of the calculation without performing a tedious 
analytic analysis. 

Some topics cannot be taught at all without computa- 
tion. The most obvious example of this is chaos. Com- 
putational solutions of systems like the driven, damped 
pendulum can be used to illustrate important concepts of 
dissipative chaos such as sensitive dependence on initial 
conditions, period doubling bifurcations, Poincare sec- 
tions, and strange attractors. Simple iterated function 
systems like the logistic map can be used to delve deeper 
into dissipative chaos with bifurcation diagrams and Lya- 
punov exponents. While these topics are now included 



FIG. 2: Trajectory of a long-range projectile fired due East 
from Rome, Georgia. One path is calculated by taking non- 
inertial forces into account, while the other path ignores these 
apparent forces. 
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FIG. 3: An OSP application for visualizing the motion of a 
charged particle in the presence of three dimensional constant 
electric and magnetic fields. 



in many textbooks) 12 * 13 ^ they cannot really be taught ef- 
fectively without letting students carry out some com- 
putations. Once students have mastered the necessary 
computational skills, even the more advanced topic of 
Hamiltonian chaos becomes accessible through the study 
of two-dimensional area-preserving maps P^j In fact, stu- 
dents who have already seen the Euler-Cromer algorithm 
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will be familiar with some of the properties of area- 
preserving maps. 

Computational projects can also be chosen to intro- 
duce important numerical methods. Numerical integra- 
tion can be introduced in the context of finding the period 
of a pendulum undergoing large amplitude oscillations. 
Numerical root-finding techniques can be used to deter- 
mine the time of flight for a projectile with air resistance. 
Numerical determination of eigenvalues and eigenvectors 
can be used in the context of finding principle axes for 
rigid-body motion or normal modes of coupled oscillator 
systems. 

Many textbooks now have computational exercises 
among their end-of-the-chapter homework problems, and 
these can be a valuable source of ideas for computational 
projects.^ One of us (JH) has written a textbook for 
classical mechanics that explicitly incorporates computa- 
tion throughout the book and includes a wide variety of 
computational projectsP^In addition, anyone can obtain 
most of our computational materials (except solutions to 
student projects) written in MatlatP^or Mathematical 
We plan to update these websites with new materials as 
they are created. We also plan to port these materials to 
Java/OSP or EJS. These materials will be made available 
through the BQ-OSP database.^ Instructors wishing to 
include computation in their physics courses are also en- 
couraged to consult computational physics texts^and the 
Open Source Physics website P 

V. SUMMARY 

A course in intermediate classical mechanics provides 
an excellent opportunity for including computation with- 
out requiring the resources and commitment needed 
for developing a new computational course or program. 
Many topics that are typically covered in an intermediate 
classical mechanics course can benefit from a computa- 



tional approach. The examples discussed in this article 
provide just a small sample. The intermediate classical 
mechanics course is also a good place to introduce compu- 
tation because it typically comes early in the curriculum. 
Ideally this course would be preceded by a computational 
physics course and followed by other upper-level physics 
courses that include computational work. We have found 
that even if computation is required only in classical me- 
chanics students will use the computational skills they 
have gained to solve problems in a wide variety of other 
physics and mathematics courses. An early exposure to 
computational work can also lead to opportunities for 
undergraduate research in computational physics. 

Each instructor must decide how much computation to 
include, which computational topics to cover, and how 
computation will be used in the course. These decisions 
must be based on a wide variety of factors, including 
the computational skills of the students, the computa- 
tional background of the instructor, the computational 
resources (both hardware and software) available, and 
time and content coverage constraints. Nonetheless, we 
feel that it is both possible and important to include 
computation somewhere in the physics curriculum. We 
hope that this article provides some guidance for those 
who wish to include computation in the curriculum, but 
cannot create a dedicated course or program in compu- 
tational physics. 
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